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The SciBooNE and MiniBooNE collaborations report the results of a i/ M disappearance search in 
the Am 2 region of 0.5 — 40 eV 2 . The neutrino rate as measured by the SciBooNE tracking detectors 
is used to constrain the rate at the MiniBooNE Cherenkov detector in the first joint analysis of data 
from both collaborations. Two separate analyses of the combined data samples set 90% confidence 
level (CL) limits on disappearance in the 0.5 — 40 eV 2 Am 2 region, with an improvement over 
previous experimental constraints between 10 and 30 eV 2 . 
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I. INTRODUCTION 

Observations of neutrino oscillations at the mass split- 
tings (Am 2 ) of ~ 10~ 5 eV 2 and ~ 10~ 3 eV 2 are consis- 
tent with three generations of neutrinos and a unitary 
mixing matrix However, the LSND 0] and Mini- 
BooNE [H experiments observe an excess of v e events 
from a beam, indicating possible new physics with 
Am 2 ~ 1 eV 2 . This could be explained by introducing 
additional generations of light neutrinos, and electroweak 
data Q require such additional light neutrino species to 
be sterile, i.e. with negligible couplings to and Z 
bosons. 

Sensitive searches for the disappearance of (and F M ) 
into sterile neutrino species can constrain models used 
to explain the LSND and MiniBooNE V e appearance re- 
sults. First, neutrino and antineutrino oscillations are 
expected to be either identical in such models or related 
via CP-violating phases Q or non-standard matter-like 
effects H, 01- Second, tests of and disappearance 

(vp— >v x ) probe elements in the neutrino mixing matrix 

that also govern v e appearance ( v^— >u e ). As a re- 
sult, a number of global analyses have been performed to 
study whether neutrino and antineutrino data sets rel- 
evant to the Am 2 ~ 1 eV 2 oscillation region in three 
different channels (v e appearance, disappearance, v e 
disappearance) can all be accommodated within sterile 
neutrino models 

The MiniBooNE collaboration has previously reported 
limits on and disappearance in this region of 
Am 2 . Substantial neutrino flux and interaction cross 
section uncertainties limited the sensitivity of the Mini- 
BooNE only analyses. Data from the SciBooNE exper- 
iment can be used to reduce these errors because Sci- 
BooNE shares the same neutrino flux and target mate- 
rial. This paper describes the disappearance analysis 
using data from both experiments. 



II. EXPERIMENTAL APPARATUS 

The muon neutrino beam, with mean energy of 
0.8 GeV, interacts in the SciBooNE (MiniBooNE) de- 
tectors at 100 m (541 m) from the neutrino production 
target. Neutrino disappearance at these distances would 
be observable as a distortion in the neutrino energy spec- 
trum and as a deficit of the total event rate. At Am 2 of 
~ leV 2 , the MiniBooNE neutrino spectrum will have a 
deficit relative to SciBooNE, where the neutrinos have 
not yet oscillated. Above ~ 1 eV 2 , both the SciBooNE 
and MiniBooNE neutrino energy spectra will be altered 
due to disappearance. 

The Fermilab Booster Neutrino Beamlinc (BNB) pro- 
vides the flux to both experiments. Protons at kinetic 



energy of 8 GeV interact with a 1.7 interaction length 
beryllium target. The mesons are then focused by a mag- 
netic field and decay in a 50 m long tunnel to produce 
the 93.6% pure v„ beam (5.9% 0.5% v e ). The BNB 
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neutrino flux [12| is simulated using GEANT4 [13j and 
includes updated p-Be particle production data jlj, [l5[ . 

The SciBooNE experimental hall houses three sub- 
detectors; in order from upstream to downstream they 
are a fully active and finely segmented scintillator tracker 
(SciBar), an electromagnetic calorimeter (EC), and a 
muon range detector (MRD). The SciBar detector [l|| 
consists of 14,336 extruded polystyrene (CsHg) strips. 
The scintillators are arranged vertically and horizontally 
to construct a 3 x 3 x 1.7 m 3 volume with a total mass of 
15 tons which serves as the primary neutrino interaction 
target for SciBooNE. Each scintillator strip is read out 
by a 64-channel multi-anode photomultiplier tube (PMT) 
via wavelength shifting fibers. The EC is a "spaghetti" - 
type calorimeter; 64 modules made of 1 mm scintillating 
fibers embedded in lead foil are bundled in 64 modules 
and read out at both ends by PMTs. The MRD is built 
from 12 iron plates, each 5 cm thick, sandwiched be- 
tween planes of 6 mm thick scintillation counters; there 
are 13 alternating horizontal and vertical planes read out 
via 362 individual 2 inch PMTs. A GEANT4 Q3] Monte 
Carlo (MC) simulation is used to model particle and light 
propagation through SciBooNE and the surrounding ma- 
terial. 

The MiniBooNE detector [l7| is a mineral oil 
Chcrcnkov light detector. It is a 12 m diameter spheri- 
cal tank lined with 1280 inward-facing 8 inch PMTs with 
an optically isolated outer region used to reject cosmic- 
ray induced events. Neutrino interactions in the oil pro- 
duce charged particles which cause emission of primar- 
ily Cherenkov light and a smaller amount of scintillation 
light. A GEANT3 [H[ based simulation, tuned on Mini- 
BooNE and external data, is used to model light and 
particle production and propagation within the detector. 

Neutrino interactions are simulated with the v3 NU- 
ANCE event generator [l9j . with cross section parame- 
ters set to match MiniBooNE and SciBooNE data 12011 . 



III. EVENT SELECTION AND 
RECONSTRUCTION 

Prior to selection, approximately 44% (39%) of all 
events in SciBooNE (MiniBooNE) are due to charged 
current quasi-elastic (CCQE) scattering and 25% (29%) 
to charged current single pion production (CCItt); the 
small difference in event fractions in MiniBooNE and 
SciBooNE is due to the somewhat higher energy flux at 
MiniBooNE. While all CC and NC interactions are 
sensitive to oscillation into sterile states, CCQE is the 
dominant interaction mode after selection cuts. The re- 
constructed neutrino energy (E® E ) is calculated assum- 
ing the neutrino interaction was CCQE, and assuming 
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the target nucleoli is at rest: 

e qe = 2(M n - E B )E ll - (E% - 2M n E B + AM + M%) 
2[(M n -E B )-E fl +p fl cose^} 

(1) 

where AM = M% — M£, M indicates the muon, proton 
or neutron mass with appropriate subscripts, E B is the 
nucleon binding energy which is set to 34 MeV, E^(p^) 
is the reconstructed muon energy (momentum), and 9^ 
is the reconstructed muon scattering angle with respect 
to the neutrino beam direction. 

Charged current (CC) neutrino interactions in 
SciBooNE are selected by identifying single muon 
tracks [2(| . The p^ and 6^ of each muon track are recon- 
structed from the hits and energy loss in the SciBooNE 
sub-dctectors. The highest momentum track in the beam 
on-time window is required to have p^ > 0.25 GeV/c 
to reduce the number of neutral current (NC) events. 
The energy loss of the track in SciBar must be consis- 
tent with a muon hypothesis, and must originate within 
the 10.6 ton SciBar fiducial volume (FV). Muon candi- 
date tracks, where the endpoint of the track is also in the 
SciBar FV, are classified as SciBar-stopped tracks; tracks 
which stop in the MRD constitute the MRD-stopped 
sample. Tracks which exit the last layer of the MRD 
form the MRD-penetrated sample. 

Figure Q] shows the event selection efficiency of each 
SciBooNE sample for CC interactions which occur within 
the SciBar FV, as estimated by MC simulation. We es- 
timate the selection cuts for the MRD-stopped sample 
are 18% efficient for all CC interactions, with 90% CC 
neutrino interaction purity; the main sources of ineffi- 
ciency stem from the detector volume acceptance and the 
momentum requirement that the events be contained in 
SciBar or the MRD, or exit through the last plane of 
the MRD. The MRD-penetrated sample selection cuts 
are 4% efficient and produce a 97% pure CC sample. 
The SciBar-stopped cuts are 12% efficient and the sam- 
ple is 85% pure CC interactions. The mean neutrino 
energies of the selected events are 1.2, 2.4 and 1.0 GeV 
for MRD-stopped, MRD-penetrated and SciBar-stopped 
samples, respectively. Together, the SciBooNE samples 
cover the muon kinematic region (p^-Op) relevant for 
events in MiniBooNE 0. 

The selection criteria yield 13589 data events in the 
SciBar-stopped sample, 20236 data events in the MRD- 
stopped sample, and 3544 events in the MRD-penetrated 
sample after subtracting cosmic ray events, for 0.99 x 
10 20 protons on target (POT) collected; with 0.3 GeV < 
E® E < 1.9 GeV, there are 13592 (20166) data events in 
SciBar-stopped (MRD-stopped) samples. 

The E® E resolution is 13% for CCQE events in 
the MRD-stopped sample, and 20% for SciBar-stopped 
events; the MRD-penetrated sample can only provide the 
muon angle for each event, not the neutrino energy. Ac- 
cording to the simulation, the SciBar-stopped sample is 
51% CCQE and 31% CCItt. The remaining events are 
from CC or NC multi-pion events in SciBar. The MRD- 
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FIG. 1. (Top) Predicted number of CC events in the SciBar 
FV as a function of true E v . The number of selected events 
in each sub-sample are also shown. (Bottom) Detection ef- 
ficiency as a function of true neutrino energy for each sub- 
sample [2(3 ]. 



stopped sample is 52% CCQE and 34% CCItt; MRD- 
stopped sample background events also include interac- 
tions in the EC or MRD which send a muon backwards 
into SciBar. The MRD-penetrated sample is 57% CCQE 
and 32% CC17T. The small background due to cosmic 
events is estimated from beam-off time windows. Fig. [5] 
shows the E® E spectrum for SciBar-stopped and MRD- 
stopped data and prediction, assuming no oscillations, 
with cosmic background subtracted. Both samples are 
divided into 16 bins for 0.3 GeV < E® E < 1.9 GeV, 
0.1 GeV wide. 

The MiniBooNE event selection is similar to the Sci- 
BooNE selection in that it is based on identifying a single 
muon in the detector . Clusters of PMT hits in time 
are categorized as sub-events within a neutrino interac- 
tion event. The timing and charge response of the PMTs 
is then used to reconstruct the position, kinetic energy, 
and direction vector of the primary particle within each 
sub-event. The muon produced from CC neutrino inter- 
actions is reconstructed as the first of two sub-events, 
with the decay electron as the second. Exactly two sub- 
events arc required and both sub-events must have mini- 
mal activity in the outer veto region. The first sub-event 
must be within the beam-on time window and have more 
than 200 inner tank PMT hits (p M > 0.25 GeV/c) to elim- 
inate electrons from stopped cosmic ray muon decays. 
The mean emission point of Cherenkov light of the first 
sub-event must be within the 442 ton fiducial volume. 
The second sub-event must have fewer than 200 inner 
PMT hits and the reconstructed vertex must be within 
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FIG. 2. E® E distribution for SciBooNE data (black) with 
statistical errors, and prediction assuming no oscillations (rel- 
atively normalized by the total number of SciBar-stopped and 
MRD-stopped events) for the SciBar-stopped sample (top) 
and MRD-stopped sample (bottom). Cosmic background is 
subtracted. Attached to the prediction are the size of the 
systematic uncertainty (shaded boxes). The predicted non- 
CCQE events (dash) events are also shown. 
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FIG. 3. Top: E$ E distribution for MiniBooNE data (black) 
with statistical errors, and prediction assuming no oscillations 
for the simultaneous fit analysis (relatively normalized by the 
total number of SciBar-stopped and MRD-stopped events). 
Attached to the prediction are the size of the systematic un- 
certainty (shaded boxes). The predicted non-CCQE events 
(dash) events are also shown. Bottom: Same as above, for 
the spectrum fit analysis, where the prediction is scaled by 
the SciBooNE corrections. 



100 cm of the endpoint of the first sub-event to be con- 
sistent with an electron from the (first sub-event) muon 
decay. Requiring one decay electron eliminates most NC 
interactions and allowing only one decay electron removes 
most CC charged pion interactions. 

MiniBooNE took neutrino events in two periods. The 
first period (Run I) comprises 5.58 x 10 20 POT, with 
190,454 events with GeV < E$ E < 1.9 GeV. The sec- 
ond period (Run II) was taken during SciBooNE's run, 
comprising 0.83 x 10 20 POT and 29269 events for the 
same energy range. The two data sets are consistent 
with each other within the time dependent uncertainties, 
such as POT normalization, but are treated separately in 
the analysis, as these uncertainties cancel between Mini- 
BooNE and SciBooNE for Run II but not for Run I. The 
MiniBooNE selection is 35% efficient with 74% purity for 
CCQE events; the resultant sample has average E® E res- 
olution of 11% for CCQE events [2l|. The rest of the 
sample is dominated by CC17T events (~75%) where the 
outgoing pion is unobserved (due, e.g., to absorption in 
the nucleus). Fig. [3] shows the the E® E spectrum for 
MiniBooNE (Run I) with the prediction, assuming no 
oscillations, in 16 bins, 0.1 GeV wide except for the first 
bin (0-0.4 GeV). 



IV. SEARCH FOR MUON NEUTRINO 
DISAPPEARANCE 

A. Analysis strategy 

Two different disappearance fits were performed on 
MiniBooNE and SciBooNE candidates. The first 
method, called the simultaneous fit, is where Mini- 
BooNE, SciBar-stopped and MRD-stopped data are fit 
simultaneously to test for oscillations in both detectors 
directly. The second method is the spectrum fit analy- 
sis, where SciBooNE data are used to extract a CC in- 
clusive rate which is applied to MiniBooNE using scale 
factors as a function of true neutrino energy. The scale 
factors are extracted from the SciBar-stopped, MRD- 
stopped and MRD-penetrated samples [2(|. The cor- 
rected MiniBooNE prediction, shown in Fig. [31 is then 
compared to data with reduced uncertainties and fitted 
for disappearance accounting for oscillation in both Sci- 
BooNE and MiniBooNE. The simultaneous fit method, 
by fitting both SciBooNE and MiniBooNE samples simul- 
taneously, is advantageous by taking all correlations into 
account between the two experiments. The spectrum fit, 
by extracting the CC interaction rate into scale factors 
which depend only on neutrino energy, loses some corre- 
lation between SciBooNE and MiniBooNE which reduces 



5 



TABLE I. CC inclusive interaction rate normalization factors, with extracted uncertainties, obtained from the SciBooNE 
spectrum fit [20l |. 



Energy region (GeV) 


0.25 - 0.50 


0.50 - 0.75 


0.75 - 1.00 


1.00 - 1.25 


1.25 - 1.75 


1.75 - 


v,j. CC rate normalization factor (fi) 


1.65 ±0.22 


1.31 ±0.11 


1.36 ±0.08 


1.38 ±0.09 


1.36 ±0.12 


0.90 ±0.09 



the power of the constraint. However, the spectrum fit 
demonstrates the reduction of the individual systematics 
at MiniBooNE due to the SciBooNE constraint, because 
it is done prior to the oscillation fit. Both techniques are 
of interest to future experiments and so are presented 
here. 



B. Oscillation prediction 

For both fits, the oscillation prediction assumes two- 
flavor — > v x disappearance characterized by one 
large mass splitting (Am 2 = Amf lk ) between the ac- 
tive, light neutrino mass states k, which participate in 
standard three-neutrino oscillation and a fourth, heavier 
neutrino state h, and one oscillation amplitude (sin 2 29 = 
4|E/ M | 2 (l r |t/ M | 2 )). 

The oscillation probability, P(v^ — > v x ), is given as 



P(v» -> v x) = sm2 26* sin 2 1.27Am 



E 



(2) 



where L[km] is the distance traveled and E[GeV] is the 
neutrino energy; Am 2 [eV 2 ] is the mass splitting and the 
mixing angle, 9, is dimensionless. The mean flight dis- 
tance of muon neutrinos detected in SciBooNE is ^76 m 
and in MiniBooNE is ^520 m. Therefore, the oscillation 
probability can be non-zero for events in both SciBooNE 
and MiniBooNE in the oscillation search region, where 
neutrino energy is ~1 GeV and the mass splitting (Am 2 ) 
is ~1 cV 2 . The predicted number of events in each re- 
constructed E® E bin in each sample, Pi(sin 2 29, Am 2 ), is 
calculated assuming all and events, both CC and 
NC, can oscillate according to a particular Am 2 -sin 2 29 
point and the true neutrino energy. The incident neu- 
trino's energy and distance from the point of its parent 
decay to detection are used from the simulation to prop- 
erly represent the range of E„, L of the SciBooNE and 
MiniBooNE samples. 

The simultaneous fit prediction, pf lm , has 48 bins 
corresponding to the SciBooNE and MiniBooNE recon- 
structed E® E bins (16 SciBar-stoppcd bins, 16 MRD- 
stopped bins, and 16 MiniBooNE bins corresponding to 
Run I or Run II only) . The prediction in each bin is nor- 
malized by the sum of the observed number of SciBar- 
stoppcd and MRD-stopped events, dk, to prediction, rife, 
assuming oscillation at Am 2 , sin 2 29: 



PI 



m (sin 2 29, Am 2 ) 



E 



32 , 



Ef=i«fc(sin 2 29, Am 2 ) 



rii (sin 2 29, Am 2 



(3) 



where k and i are indices over reconstructed E® E bins. 

For the spectrum fit, the prediction p s i pec corresponds 
to 32 MiniBooNE reconstructed E® E bins (16 Run I bins 
and 16 Run II bins). The prediction in each bin is built 
by additionally applying the extracted rate normaliza- 
tion factors, fi, obtained by the a naly sis of charged 
current interactions at SciBooNE [2fJ. The normaliza- 
tion factors, fi, given in Table U as a function of incident 
neutrino energy, will be referred to as the SciBooNE cor- 
rections. As with the simultaneous fit, the prediction is 
further scaled by the expected number of oscillated events 
at SciBooNE: 



pf ec (sin 2 29, Am 2 ) 
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E 



fi 



Mj(0,0) 



A/} (sin 2 29, Am 2 ) 



Mij{sm 2 29,Am 2 ), (4) 



where A/} is the predicted number of events in true E v 
bin j at SciBooNE, and A4ij is the predicted number of 
events at MiniBooNE for reconstructed E® E bin i and 
incident true neutrino energy bin j, assuming oscillation 
with parameters Am 2 , sin 2 29. 



C. Systematic Uncertainties 

Systematic uncertainties are included for the underly- 
ing neutrino flux prediction, neutrino interaction cross 
section, and detector response. The method used to esti- 
mate the uncertainties due to the underlying neutrino 
flux prediction and the MiniBooNE detector model is 
identical to the method used in previous MiniBooNE 
results [ll], HH ; the cross section uncertainties and Sci- 
BooNE detector uncertainties assumed are identical to 
previous SciBooNE results [2(| • The dominant uncertain- 
ties on the neutrino event rate at SciBooNE and Mini- 
BooNE are the uncertainty of the neutrino flux due to 
production of 7r + mesons from the Be target, and the 
CCQE and CC17T neutrino cross sections. The detector 
uncertainties affect the muon energy and angular resolu- 
tion. For SciBooNE, these include uncertainties on the 
energy loss of muons through scintillator and iron, light 
attenuation in the WLS fibers and PMT response; for 
MiniBooNE, detector uncertainties include particle and 
light propagation through the mineral oil, scattering, de- 
tection and PMT response. Table fll] shows the ratio of 
predicted to observed number of events for the SciBooNE 
and MiniBooNE samples, and the size of the flux, cross 
section and detector uncertainties. 
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TABLE II. Ratio of observed number of events to predicted for each sample with uncertainties due to statistics, neutrino flux, 
cross section and detector systematics. 



Sample 


Data/MC 






Uncertainties (%) 






Stat 


Flux 


Cross section 


Detector 


Total sys. 


SciBar-stopped 


1.25 


1.0 


6.2 


15.3 


2.2 


16.6 


MRD-stopped 


1.26 


0.8 


8.0 


16.8 


2.8 


18.3 


MRD-penetrated 


1.04 


1.7 


18.2 


17.1 


5.8 


25.7 


MiniBooNE 


1.21 


0.3 


6.2 


15.1 


4.8 


17.1 



TABLE III. Ratio of observed number of events to predicted for MiniBooNE sample with uncertainties due to statistics, 
neutrino flux, cross section and detector systematics, after applying the SciBooNE corrections listed in Tab. [I] The effect of the 
uncertainties associated with the extraction of the SciBooNE corrections (fi) are also shown. 











Uncertainties (%) 






Sample 


Data/MC 


Stat 


Flux 


Cross section Detector 


SB spectrum 


Total sys. 


MiniBooNE 


0.90 


0.2 


0.8 


4.7 4.8 


1.7 


7.3 



TABLE IV. Ratio of observed number of events to predicted 
for SciBooNE samples after re-applying the SciBooNE correc- 
tions listed in Tab. [I] 

Sample SciBar-stopped MRD-stopped MRD-penetrated 
Data/MC 0.94 0.97 0.98 



Each source of systematic uncertainty produces corre- 
lations between E$ E bins within a single data sample 
and between samples, e.g. between SciBooNE and Mini- 
BooNE. These correlations are represented by developing 
a covariancc matrix, Vij , which is the sum of the individ- 
ual covariance matrices for each source of systematic un- 
certainty. Each covariance matrix is constructed in the 
same manner as in previous SciBooNE CC inclusive anal- 
ysis [2(| and MiniBooNE oscillation analyses [13, HH ; a 
fractional error matrix is generated using the MC expec- 
tations for the central value (pi(j)) and many different 
sets for systematic variations (p*^ k ) for each E® E bin: 

M sys sys 

v --m\— J— (5) 

where k denotes the index for each systematic varia- 
tion and M denotes the total number of variations. For 
the simultaneous fit the covariance matrix ranges over 
E® E bins from SciBar-stopped, MRD-stopped and Mini- 
BooNE samples. The sizes of the diagonal elements of 
the covariance matrix, Oi = Pi^/Vu are shown as the 
error bands on the predictions for the individual sam- 
ples in Figs, rj] and [3J and also in Figs. [7] and [FJ The 
correlation matrix component of the covariance matrix, 
/ y/VuVjj , is shown in Fig. [4] Where there is a shared 
source of uncertainty (e.g., neutrino flux or cross section), 
MiniBooNE bins are strongly correlated with SciBooNE 
bins. 



The correlations between SciBooNE and MiniBooNE 
are exploited by the simultaneous fit y 2 , described in 
Sec. IIVD1 to constrain the rate at MiniBooNE. When a 
bin at MiniBooNE is strongly correlated to one at Sci- 
BooNE, a difference between the data and prediction at 
MiniBooNE will only increase the \ 2 if the- difference is 
not also present in the appropriate SciBooNE bin. The 
differences in the event rates in the two detectors due to 
flux or cross section systematics, unlike detector system- 
atics, tend to cancel. As a result, the largest uncertainty 
for the simultaneous fit analysis is the MiniBooNE de- 
tector uncertainty. 

In the spectrum fit method, uncertainties of neutrino 
flux and cross sections are constrained by the SciBooNE 
corrections. The systematic variations of the predicted 
number of events in MiniBooNE after the SciBooNE cor- 
rections, p s .P ec > s y s ^ are calculated as: 

3 3 

where Afj and A/J ys are the central value and the system- 
atic variation (sys), respectively, of the predicted num- 
ber of events in SciBooNE, and Mlj S is the system- 
atic variation of the predicted number of events at Mini- 
BooNE for E® E bin i and incident neutrino energy bin 
j. Then, p s P ec < sys i s used to construct the covariance ma- 
trix for the spectrum fit analysis for MiniBooNE E® E 
bins. The ratio of the observed number of events pre- 
dicted for the MiniBooNE sample and the size of the 
reduced systematic uncertainties after applying the Sci- 
BooNE corrections are shown in Tab. IIIII As a closure 
test, the SciBooNE corrections are re-applied to the Sci- 
BooNE MC samples and compared with the data, as 
shown in Tab. IIVI The data to MC ratios for the Sci- 
BooNE samples are not exactly zero due to the small 
data and MC disagreement at the low Q 2 region [2fj| . 
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Reconstructed E v bin number Reconstructed E^GeV) 



FIG. 4. Correlation coefficients of the total systematic un- 
certainties on the reconstructed E v distribution for the simul- 
taneous fit analysis. Bins 0-15 and bins 16-31 are respectively 
for the SciBar-stopped and the MRD-stopped samples from 
SciBooNE, shown in Fig.d Bins 32-47 are for the MiniBooNE 
sample shown in the top panel of Fig. 



FIG. 5. The fractional size of the systematic uncertainties 
for each MiniBooNE reconstructed E v bins for the spectrum 
fit analysis. The dashed lines show the uncertainties using 
MiniBooNE alone, and the solid lines show the uncertainties 
after the SciBooNE corrections. 



however the differences are well within the systematic 
uncertainties. The fractional sizes of the diagonal ele- 
ments of the covariance matrix are shown in Fig. [5J and 
the correlation coefficients for the total systematic un- 
certainty after the SciBooNE corrections are shown in 
Fig. [6] for MiniBooNE E$ E bins. After the SciBooNE 
corrections, the flux and cross section uncertainties are 
reduced by a factor of approximately two as compared to 
the previous MiniBooNE-only analysis. The spectrum fit 
method has uncertainties associated with the extraction 
of the fi rate factors, given in Table U which are included 
in Vij. As shown in Tab. Hill the dominant uncertainty 
for the spectrum fit analysis is the MiniBooNE detector 
uncertainty which does not cancel in the ratio between 
the two experiments. 

D. Oscillation Fit 

The disappearance search uses a frequentist-based 
A% 2 approach [24| to determine allowed regions in the 
Am 2 — sin 2 29 plane. The \ 2 is calculated from a com- 
parison of the data, di, in E® E bin i, to a prediction 
K (Am 2 , sin 2 29): 

rr* Pi Pj 

where Vij is the covariance matrix. As defined in Sec- 
tion IIVB1 the simultaneous fit prediction pf im has 48 




Reconstructed E v (GeV) 

FIG. 6. Correlation coefficients of the total systematic uncer- 
tainties on the MiniBooNE reconstructed E v distribution for 
the spectrum fit analysis, shown in the right panel of Fig. [3] 



bins and the spectrum fit prediction, p| pec has 32 bins. 
The covariance matrix Vij is either 48 x 48 bins (simul- 
taneous fit) or 32 x 32 bins (spectrum fit). At each 
oscillation point in a grid across the region of interest, 
Ax 2 (Am 2 , sin 2 29) = \ 2 (Am 2 , sin 2 29) ~x 2 ( min ) is cal- 
culated. The probability distribution for the A% 2 statis- 
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0.4 0.6 0.8 1 1.21.41.61.8 0.40.60.8 1 1.21.41.61.8 0.20.40.60.8 1 1.21.41.61.8 

Reconstructed E, (GeV) Reconstructed (GeV) Reconstructed (GeV) 



FIG. 7. Ratios with respect to the null hypothesis for SciBar-stopped (left), MRD-stopped (middle) and MiniBooNE (right) 
samples for the simultaneous fit analysis. Data (points), and disappearance expectations for the best fit (solid red) and 
parameters Am 2 = 1.0 eV 2 , sin 2 28 = 0.5 (dashed blue) and Am 2 = 10.0 eV 2 , sin 2 26 = 0.5 (dotted magenta) are shown. The 
best fit parameters are Am 2 = 43.7 eV 2 , sin 2 26 — 0.60. Fractional size of the diagonal elements of the error matrix for the 
null oscillation is shown as the gray bands. 



S. 1-4 
1.2 

1 

0.8 
0.6 



Best fit 

Am 2 =1.0eV 2 , sin 2 20=O.5 
Am 2 =10.0eV 2 , sin 2 20=O.5 

+ 




0.20.40.60.8 1 1.21.41.61.8 
Reconstructed E v (GeV) 



FIG. 8. Ratios with respect to the null hypothesis for Mini- 
BooNE samples for the spectrum fit analysis. Data (points), 
and disappearance expectations for the best fit (solid red) and 
parameters Am 2 = 1.0 eV 2 , sin 2 26 = 0.5 (dashed blue) and 
Am 2 = 10.0 eV 2 , sin 2 26 = 0.5 (dotted magenta) are shown. 
The expectations are for the spectrum fit analysis. The best 
fit parameters are Am 2 = 41.7 eV 2 , sin 2 26 — 0.51. Frac- 
tional size of the diagonal elements of the error matrix for the 
null oscillation is shown as the gray bands. 



tic is determined from an analysis of a set of simulated 
data samples generated with neutrino flux, cross section, 
and detector systematics varied. 



V. RESULTS AND CONCLUSION 

Both the simultaneous fit and spectrum fit are consis- 
tent with the null hypothesis. Fig. [7] and Fig. [3] show 
the ratio of data to the null hypothesis and three oscilla- 
tion scenarios for illustration, for the simultaneous fit and 
the spectrum fit, respectively. The x 2 value for the null 
hypothesis is 45.1 (48 degrees of freedom (DOF), 59% 
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90% CL limits from CCFR and CDHS 
90% CL limit from MINOS 
MiniBooNE only 90% CL sensitivity 
MiniBooNE only 90% CL limit 
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90% CL observed (Sim. fit) 
0% CL observed (Spec, fit) 
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0.9 1 
sin 2 2 6 



FIG. 9. The 90% CL limit for the simultaneous fit (solid 
black curve), and for the spectrum fit (red dashed curve). The 
sensitivity for the simultaneous fit analysis is also shown (dot- 
dash curve) with previous limits by CCFR [2f| and CDHS I2I 
(light gray), MiniBooNE [TJ] (dark gray) and MINOS O 
(green hash box). The sensitivity for the MiniBooNE-only 
analysis [XiTJ is also shown (cyan dashed curve). 



probability), with x 2 ( m i n ) = 39.5 at Am 2 = 43.7eV 2 , 
sin 2 29 = 0.60 in the simultaneous fit analysis, fit to Mini- 
BooNE Run I data. The expected Ax 2 for the 90% CL 
for the simultaneous fit is 9.34, so the result is consistent 
with the null hypothesis. A fit to Run II data is consistent 
with a X 2 (null) = 41.5 (48 DOF, 73% probability); a fit 
to both Run I and Run II produces a negligible improve- 
ment on the simultaneous fit sensitivity as compared to 
fitting just the larger Run I data set. For the spectrum 
fit to Run I and Run II, the x 2 found comparing the data 
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to the null hypothesis is 41.5 (32 DOF, 12% probability) 
with the x 2 minimum value of 35.6 at Am 2 = 41.7eV 2 , 
sin 2 26> = 0.51. The expected A X 2 for the 90% CL is 8.41, 
so the fit is consistent with the null hypothesis. 

Figure |9] shows the 90% CL sensitivity for the simulta- 
neous fit analysis compared to the previous MiniBooNE- 
only sensitivity. The 90% CL sensitivity is the average 
limit of 1000 fake experiments with no underlying oscil- 
lation assumed. Both sensitivities use the MiniBooNE 
Run I data to show how the addition of SciBooNE data 
improves the sensitivity. In the MiniBooNE-only fit, flux 
and cross section uncertainties are dominant at low Am 2 
(Am 2 < leV 2 ), resulting in a reduced sensitivity in that 
region. In the simultaneous fit, the presence of flux and 
cross sections uncertainties does not increase the x 2 once 
SciBooNE data arc included and the sensitivity improves. 
In the spectrum fit, which has a similar sensitivity as the 
simultaneous fit, the uncertainties are directly reduced 
by the SciBooNE corrections and so the sensitivity is im- 
proved. 

The 90% CL limit for both analyses is also shown in 
Fig. [S] along with other relevant experimental results. 
While the limit is an improvement in the Am 2 region 
of 10 — 30 cV 2 , the 90% CL limits are worse than the 
sensitivity at Am 2 ~ 1 eV 2 because of the small deficit 
at MiniBooNE relative to SciBooNE. The energy depen- 
dence of this deficit is within the uncertainties, thus the 



limit is consistent with the region covered by 68% of the 
limits of fake null experiments. The spectrum fit limit 
is not as strong as the simultaneous fit limit because the 
corrected spectrum creates a larger deficit between Mini- 
BooNE data and prediction than the default prediction, 
while the simultaneous fit applies no correction and so is 
consistent with the small deficit. 

In summary, we search for disappearance using both 
SciBooNE and MiniBooNE data sets. The two analyses 
presented reduce the flux and cross section uncertainties 
to the level of the detector uncertainties of the two exper- 
iments. We set 90% CL limits on disappearance in the 
Am 2 region of 0.5 - 40 eV 2 , with an improvement over 
previous experimental results between 10 and 30 eV 2 . 
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